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Abstract 



We study the normalized trace g n (z ) = n _1 tr (H — zl)^ 1 of the resolvent of n x n 
real symmetric matrices H = [(l + 5jk)Wjk/\/n\j k=1 assuming that their entries 
are independent but not necessarily identically distributed random variables. 
We develop a rigorous method of asymptotic analysis of moments of g n {z) for 
\^ z \ > Vo where r/o is determined by the second moment of Wjk- By using 
this method we find the asymptotic form of the expectation ~Ei{g n (z)} and of the 
connected correlator Ei{g n (zi)g n (z2)} — E{g n (zi)}~E{g n (z2)}. We also prove that 
the centralized trace ng n (z) — ~E{ng n (z)} has the Gaussian distribution in the 
limit n = oo. Basing on these results we present heuristic arguments supporting 
the universality property of the local eigenvalue statistics for this class of random 
matrix ensembles. 

PACS number (s): 05.40+j 
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I. INTRODUCTION 



Since the pioneer works of Wigner and Dyson random matrix theory (RMT) has been 
successfully used to describe the energy levels of complex quantum systems: heavy nuclei, 
quantum chaotic systems, mesoscopic samples, etc (see e. g. Refs. |H|). Another rather broad 
field of the RMT applications is related to quantum field theory: large colour limit of QCD, 
2D quantum gravity and bosonic strings (see e.g. Refs. 

The phenomenological nature of the RMT approach that may be regarded as its certain 
drawback on the one hand, provides, on the other hand, the model independent frameworks, 
that make the approach applicable to a wide variety of systems having different microscopic 
natures and origins. These frameworks assume certain amount of "robustness" of the RMT 
models and results. In other words, it is believed that "sufficiently large" number of them 
should have no dependence or a rather weak one on a random matrix ensemble used. This 
belief explains partly the fact that majority of the RMT ideas and applications are based on 
results obtained for the archetype Gaussian Ensembles (GE's) and the Circular Ensembles 
(CE's). On the other hand, this belief requires certain justification, in particular extending the 
results known for the GE's and for the CE's to other classes of ensembles. 

The most referred to are the Gaussian Orthogonal Ensemble (GOE) of random n x n 
symmetric matrices and the Gaussian Unitary Ensemble (GUE) of random n x n Hermitian 
matrices. The density of the probability distribution in these ensembles has the form 

P n (H) = Z- 1 ex V [-ntrF(H)}, (LI) 

where F(x) = x 2 /Aw 2 and Z n is the normalization constant. 



The probability distribution (|I.1|) possesses two important properties: i) it is invariant with 
respect to either orthogonal or unitary transformations of R n or C n , respectively; and ii) the 
matrix elements are independent random variables (modulo the obvious symmetry conditions). 

These properties of the GE's determine them uniquely and motivate two classes of general- 
isations of the GE's. 

The first class consists of ensembles having an orthogonal or unitary invariant but not 
necessarily matrix-element-independent probability distribution. The typical representatives 



are the ensembles with the probability distribution of the form (|1.1[) in which F(x) is an arbitrary 



bounded below and growing fast enough on infinity function. These invariant ensembles can 
be used to describe physical systems having no preferential basis. They arose also in studying 
the large-n limit in quantum field theoryEHj and later found other applicationslHlll. 

Random matrices with invariant distributions show remarkable "robustness" (known as the 
universality) of spectral properties in the microscopic regime. In this regime one scales the 
energy so that the mean distance between nearest eigenvalues remains of order unity as the 
dimension of matrices increases&0. Thus one is able to study properties of a finite number 
of eigenvalues. The universality of the level spacing distribution and other microscopic (lo- 
cal) spectral characteristics has been extensively discussed in recent theoretical physics and 
mathematical literature. We refer the reader to a number of publications: Refs. pTl|— 116|. 

The second class consists of ensembles whose matrix elements in a certain basis are inde- 
pendent random variables, i. e. the ensemble probability distribution factorizes into a product 
of distributions of the matrix elements in this basis. The corresponding random matrices can 
be associated with physical systems having a preferential basis and appear, in particular, in 
condensed matter physics and theory of disordered systems. This second class goes back to 
Wigne and we shall refer to the corresponding ensembles as Wigner ensembles (or Wigner 
matrices) . 

The subject of the present paper is the Wigner ensemble of n x n real symmetric matrices 
of the form 



H = [H jk \l k=1 , H jk = (1 + 5 jk )W jk /yfti, (1.2) 
where Wj k ,j < k are independent random variables such that 

E{W jk } = 0, E{Wf k } = w 2 . (1.3) 

Here and thereafter E{-} denotes averaging over all Wj k ,j < k. 

The distributions of Wj k s may depend on (j, k), but we assume that they are independent 
of n. We make the latter assumption mainly for the sake of technical simplicity. On the other 
hand, this assumption allows one to consider all Wj k on the same probability space and to find 
an optimal form of a number of important facts related to the Wigner ensembles (for example, 
the convergence with probability 1 in formulae ([L7]) and (|I.11|) below). If Wj k s are independent 
Gaussian random variables, then the ensemble (|E2])-(|rT3"D coincides with the GOE. This justifies 
the presence of the term with 5j k in ([L2j). 

Macroscopic properties of Wigner ensembles are more or less well understood. We call 
macroscopic the asymptotic regime in which the number of eigenvalues in unit energy interval 
is proportional to n. Discussing macroscopic properties of random matrices we have to mention 
first of all the density of states (DOS) which is the simplest macroscopic characteristic of the 
ensemble eigenvalue statistics. It turns out that under rather natural and mild conditions on 
the distributions of Wj k the DOS in the Wigner ensemble ([Ti2|)-(|rT3"D does not depend on the 
form of the distributions of Wj k . This DOS is known as the Wigner semi-circle law (see Eq. ( |L5| ) 
below) . Other macroscopic spectral quantities such as the conductivity and the interband light 
absorption coefficient show the same "robustness" Definition of these quantities requires 
some care for Wigner matrices. However, the conductivity and the interband light absorption 
coefficient can be defined and computed for the so-called band random matrices and random 
operators with independent matrix elements that are quite close to the Wigner ensemble ( |L2"|) - 
( |Q| ) in their macroscopic properties both technically and by results ( see e.g. Ref. f^). As 
for the microscopic scale, supersymmetry calculations^ suggest "robustness" (universality) of 
spectral properties of the Wigner ensemble (|L"2"|)-([0|) as well but evidence of this has not been 
rigorously established so far. 

Introduce the normalized eigenvalue counting function 

N n (E) = n^ftiEj : Ej is an eigenvalue of H and Ej < E}. (1.4) 

Wigner in the end of fifties proved0 that in the case of identically distributed Wj k having all 
moments N n (E) converges in probability as n —>■ oo to a non- decreasing function N SC (E) (the 
semicircle law) whose derivative (DOS) is 

fjr\ \ V4w 2 - E 2 \E\ < 2w , TK x 
P(E) = < 2ttw 2 1 1 ~ (1.5) 

I \E\> 2w. 

The modern formulation of Wigner's result is as follows. Let us consider random matrices 
([L2])- (lL3D with mutually independent arbitrary distributed entries defined on a common proba- 
bility space. Then the condition (the matrix analogue of the Lindeberg condition of probability 
theory) 



is 



lim — V / x 2 dProb\W ik < x] = 0, for any v > 0, (1.6) 

sufficientii and necessary^ for the following limiting relation 

lim N n (E) = N SC (E). (1.7) 



to hold for every E with probability If we will not assume that Wjk are defined on the 
same probability space or if their probability distributions depend on n, then the same condition 
( |L6| ) will imply the convergence in probability in ( |I.7|) 

As usual in spectral theory, this result admits a natural reformulation in terms of the 
resolvent (Green's function). Indeed, the normalized trace of the resolvent 

g n (z) = n-hr (H - zl)- 1 (1.8) 

is simply the Stieltjes transform of N n (E): 

n jt x Ej — z J E — z 



Denote the Stieltjes transform of the Wigner law ( |1.5| ) by 



r[z) 



N sc (dE) -z + yg - 4w* 



The obvious condition ^sr(z)Qz > determines the branch of the square root in (|I.10| ). Due 
to one-to-one correspondence between non-decreasing functions and their Stieltjes transformsil 



( H7 ) is equivalent to the following limiting relation 

lim g n (z) = r(z), (1.11) 

n— >oo 

which holds with probability 1 for any non-real z. 

The relation ( [1.1 1| ) and the obvious bound 1^(^)1 < l^- 2 !" 1 imply that the variance of g n (z) 
vanishes as n — > oo and hence the moments 

v 



m<p{z 1 ,..., Zp )=E{i[g n (z l )} (1.12) 

factorize: 



i=i 



m 



W(z 1 ,...,z p )=i[m<p(z l )+o(l), rwoo (1.13) 



i=i 



This factorization, that follows already from the convergence in probability in (|T77| ) or in ( L.ll ) 



is typical for the macroscopic regime and can be found hidden behind many calculations in 



this regime. It is known in fact since papers Refs. [T7|, |27|, |23] and |28. 

Since according to (11) m^'(z) = E{g n (z)} = r(z)+o(l), the leading term of m^(zi, . . . , z p ) 
is Ilf=i f(zi). It seems interesting from a number of points of view to find also sub-leading terms 
and their dependence on the probability distributions of matrix elements. For instance these 
sub-leading terms are important when we would like to go beyond the macroscopic regime, 
when we are computing connected correlators of g n (z), etc. 

For the Gaussian entries respective corrections were studied in Ref. ^9] where the formal per- 
turbation theory with respect to and the respective diagrammatic technique were applied. 
This approach is an adaptation of the technique developed in Ref. ^ in order to construct the 
1/n expansion for the random operator describing disordered systems on Z d with n orbitals per 
site. 

In Ref. |3T] we suggested an approach that allows for the rigorous treatment of this problem 
in the general case of independent and arbitrary not necessarily identically distributed matrix 
elements. Our approach allows us to estimate remainders in respective asymptotic formulae 
and we show that these estimates are in a sense optimal. The approach is also free to the 



large extent from the cumbersome combinatorial problem of rearranging diagrams which is 
neccessary in order to carry out various "dressing" procedures. In particular, the dressing 
procedure that replaces the "bare" Green function —1/z by lim^oo E{g n (z)} is automatic in 



our approach. Following Ref. |31] one is able to find as many terms in the asymptotic expansion 



of m^(zi, . . . , z p ) as needed, though the technical difficulties increase with the order. 



In the present paper we use the general scheme of Ref. |31] in order to compute first terms 
in the asymptotic expansion of E,{g n (z)}. We also prove that if the distributions of Wjk satisfy 
the Lindeberg condition ( |L6|) with x 2 in the integral replaced by x 4 , and if in addition to ( p3|) 
the fourth moments of Wjk do not depend on (j, k), then 



where 



Z-2 



m 



(2), 



Zl,Z 2 



m { V{z l )mV{z 2 )=n~ 2 f{z l ,z 2 ) + o{n- 2 ), 



(1.14) 



f{zi,z 2 ) 



(1.15) 



2w 2 



[1 — w 2 r 2 (zi)][l — w 2 r 2 (z 2 )} z\—z 2 



rizx) - r(z 2 



+ 



2ar 3 (z\)r 3 (z 2 ) 



[1 — w 2 r 2 (zi)][l — w 2 r 2 (z 2 )} 



and a = E{W 4 k } — 3E 2 {H^ 2 fc } is the excess of Wjk- We also establish a limit theorem for the 
centralized trace ng n (z) — E{ng n (z)} of the resolvent. 

Unfortunately, our approach gives a bound for the remainder term in ( |L14|) containing a 
power of IQ^Q^I -1 as a factor. Thus we cannot treat rigorously the microscopic regime which 
requires Qz oc 1/n. On the other hand, the first term (15) of the asymptotic formula (14) is 
well defined in this regime and coincides with respective exact expression known for the GOE, 
provided that the latter is considered for large level spacings and is smoothed over an interval 
A such that l/n^c|A|^Clin proper units (see Section VI). We feel therefore, that by using 
our procedure of the computing of corrections, i.e. keeping the imaginary part of energy fixed 
when n goes to infinity and then letting Qz go to zero , one may treat energy intervals that are 
very large on the microscopic scale. On this intermediate scale the second term in the r.h.s. of 
(|I.15| ) which contains the probability distribution excess o vanishes and the above mentioned 
universality restores. 

Our article is organized as follows. In section || we present our basic tools. In sections |T| we 
calculate first terms of the asymptotic expansion for E{g n (z)}. In section [TV] we give a simple 
proof of ( |I.14|) - (|I.15|) with o(n~ 2 ) replaced by 0(n~ 5 / 2 ) provided that 5th absolute moment of 
Wjk is uniformly bounded. This result was cited without proof in Ref. 51. In section |V| we 
treat the general case of Wjk satisfying the higher order Lindeberg condition mentioned above. 
We prove that the fluctuations of ng n (z) around its mean value become Gaussian in the limit 
n — > oo and that the covariance of the limiting Gaussian function is f(zi,z 2 ), thus proving 
(|I.14j )-( I.15 ) in the general case. Section VI contains a discussion of some implications of our 
results. 



II. PRELIMINARIES 

In this section we present our basic technical tools. 

(i) If £ is a real-valued random variable such that E{|£| p+2 } < oo and if f(t) is a complex- 
valued function of a real variable such that its first p + 1 derivatives are continuous and 
bounded, then 

mm)} = E E{/W(0} + e, (IL16) 

a=0 a - 



where n a are the semi-invariants (cumulants) of £, |e| < C sup t \f <yP+l \t)\ E{|£| p+2 } and 
the quantity C depends on p only. 

The semi-invariants can be expressed in terms of the moments. If E{£} = (the case we 
shall deal with) and fi a = E{£ a }, then few first such relations are: K\ — fMi — 0, n 2 = fi 2 , 
K 3 — A*3; ^4 — A*4 — 3/U,|, k 5 = /i 5 — 10/^3/x 2 , ^6 = ^6 — 15^4/i2 — lO/i 2 — 30/i2, etc. For a 
Gaussian random variable with zero mean, all semi-invariants but k 2 vanish and ([II.16|) 
reduces to the exact relation 

E{£/(0> = E{^ 2 }E{/'(0}, (11.17) 

which can directly be checked integrating the l.h.s. of (|II.17| ) by parts. This is only 
the case , when formula ( Ll.l(j| ) contains finite number of terms for non-polynomial /'s. 



Indeed, according to the Marcynkiewicz theoremU if all but finite number of cumulants 
are zero, then only first and second can be nonzero. 

(ii) For any matrix A = {A aP ]^ =1 

° U- 1 ) =-(A~ 1 ) (A-i 



dAjk ^ ' a P ^ ' aj \ ) k/3 

provided A -1 exists. For the resolvent G of a real symmetric matrix H this becomes 

(11.18) 



dG a _ J —G a jG k/ 3, j = k 



dHj k [ —G a jGkf3 — GukGjp, j k 

(iii) For any two real symmetric matrices and any non-real z the resolvent identity 

(H 2 - zl)- 1 = (#! - zl)- 1 - (Hi - ziy\H 2 - #0(F 2 - ziy 1 (11.19) 

is valid. In particular, if H 2 = H, Hi = and G = {H — zl)^ 1 , then 

n 

Gj m = z 5j m + z ^ GjkH km (11.20) 



k=l 



Let H belong to the Wigner ensemble (|1.2|)-(|OD. For a fixed complex z consider complex- 
valued random variable g n (z) = n _1 tr (H — zl)~ l . Define its variance as 

E{\g n (z) - E{g n (z)}\ 2 } = F n (z, (11.21) 

and define also the domain in the complex plane as follows 

U = {z E C± : \Qz\ > 2w}. (11.22) 

We use t to denote complex conjugate and the sub-index C to indicate centering to zero 
mean. For instance, <7„ (z) = g n { z ) — E{g n (z) and thus we can rewrite ( |11.21| ) as 



F n (z,J) = EUg^z)] 2 } = E{g n (z)g^)}. (11.23) 

We will write 0(n~ p ) in asymptotic formulae for the remainders having an uniform (with respect 
z G Uq) upper bound of the form Cn~ p where C does not depend on n. In fact the bounds 
we are able to derive contain 1/(1 — \$sz\ 2 /2w 2 ) (see e.g. formula (29) below for the simplest 
case). Thus C is finite for any fixed z satisfying \Qz\ > y2w. But we prefer to use \$sz\ > 2w 
in favour of uniformity of the bounds with respect to z £ Uq. 



(iv) Let H belong to the Wigner ensemble (|L2"|)-([L3|). Assume that the 5th absolute moment 
of the random variables Wjk is uniformly bounded, i. e. sup J<fc EjlW^-fcl 5 } < +00, and 
that z EUq. Then 



E{\g n (z) - B{g n (z)}\ 2 } = E{\g^(z)\ 2 } = 0(n~' 2 ), as n - 00 (11.24) 



Let us comment on (i)-(iv). Facts (ii) and (iii) are well known. The "decoupling" formula 



( [II.16|) is simple to understand in the case when £ has all moments and f(x) belongs to the 
Schwartz space. Indeed, by using the Parseval relation for the Fourier transforms we can rewrite 
the l.h.s. of flm6D as 



j poo A j poo A 

S L Jt F ^ dt = -SF L ^)s n W* < n ' 25 > 

where 

/oo roc 
e*7(0# and U(t) = / e^'dP® 
-00 J— 00 

are the Fourier transforms of /(£) and of the probability distribution P(-) of £, respectively. 
Now, if we take into account that 

m = g and u{t) = log n ( t) = £ 

a=0 °" a=l a ' 



we can rewrite the r.h.s. of ( |11.25| ) as 



-f r F\t)u\t)e^dt = -f E^ff (^) a+1 Ft(t)n(t)rft 

Z7T J-00 Z7T ^Tq a! J-00 

00 

= E^T i E{/W(£)} 



a=0 fl! 

where we again used the Parseval relation. The latter formula is obviously (|II.16|) for p = 00. 
The case when £ has a finite number of moments and /(£) has respective number of derivatives 
requires certain technicalities which we will not discuss here. 

The bound (iv) plays an important role in many questions of random matrix theory and 
its applications. In fact, it is the simplest of bounds for connected correlators (cumulants) 
of g n (z) or, more generally, for cumulants of linear statistics of the eigenvalues ( i.e. sums 
n l £j=o 4>{Ej) where (f>(E) is smooth enough). We are going to present a detailed derivation of 
these bounds and asymptotics (both, for the Wigner ensembles and unitary invariant ensembles) 
in a subsequent publication. 

Here we only outline the scheme of the derivation of the bound (iv) considering the simplest 
case of the GOE and treating it as a representative of the Wigner ensembles, i.e. ensembles 
with independent entries. Set r n (z) = E{g n (z)}. Then, according to ([11.171) , ( |1L18|) and (jll.20| ) 
we have the relation 

r n (z) = V{g 2 n (z)} + — E{tr G 2 } (11.26) 

z z n z z 



Applying similar arguments to E{|^(2;)|} and using ( [11.261 ) we obtain the analogous relation 
for the variance ([11.23| ) 



= --E{g 2 n {z)g c n ^)} - — E{[^(z)]ttr G 2 } - — E{tr G(G*) 2 } (11.27) 



. o L L-?n V /J J q 

2 77/2: n d z 



where G* = (H - z^I)" 1 . By using the identity E{g 2 n (z)g^(z^)} = E{(g n (z) + 
E{g n (z)})\g^(z)\ 2 }, Cauchy-Schwarz inequality and the inequality ( |III.33|) below we can show 
that the first term in the r.h.s of ( |11.27p is bounded above by 2w 2 r]~ 2 F n , where rj = \5sz\, the 
second is bounded by w 2 (ni] 2 )~ 1 F^ 2 and the third is bounded by 2w 2 (n 2 ?7 4 )~ 1 . As the result 
we obtain the following inequality for rj 2 > 2w 2 



2w 2 . ^ w 2 * /9 2w 2 

-)F n - —F X J 2 - — A < 0. (11.28) 



which implies that 



J Ti o n oa 

7] z nr] z n z rf 



E{|^(z)| 2 } = E{\g n (z) - E{g n (z)}\ 2 } < % (11.29) 



where G\ = ?7~ 2 (e - l)~ 2 Ci(e), e = r] 2 /2w 2 , 1 < e < oo and Ci(e) is finite for 1 < e < oo 
Thus we have obtained ( |11.24| ) for the GOE. This is a simplest but typical bound that can 
be obtained by our method. In the general case of non-Gaussian Wjk's one has to iterate the 
resolvent identity (20) and use (16) instead of (17), truncating this procedure on the proper 
step and estimating remainders by variants of arguments presented above. 

The bound ( |1T29| ) and ( PX26D allows us to prove Q and (pT|) for the GOE. Indeed, 
combining (|II.26|) and ( |IL29|) we obtain 

/ s 1 W 2 / m Go 

where C2 has same properties as C\ in ( |11.29| ). This bound and standard compactness arguments 
show that any limit point r{z) of the sequence {r n (z)} satisfies the equation 

w 2 r 2 (z) + zr{z) + 1 = 0. (11.30) 

for \^sz\ > i] > 0. Since this equation has the unique solution ( |I.10j ) satisfying ?sr(z)^sz > 0, we 
conclude that uniformly in \5sz\ > i] > lim^oo r n (z) = r(z) where r(z) is given by (|I.10|) . 



Besides, since the Gaussian Wjk satisfying (|IT3"D can be defined on the same probability space 
we conclude from ([II. 29 ) and the Borel-Cantelli lemma that (|I.11|) and ( p7|) and are valid. 



III. ASYMPTOTIC EXPANSION FOR E{g n (z)} 

We recall our notation m^{z) for the mean value of g n {z) = n _1 tr (H — zl)" 1 . In this 
section we prove the following 

Theorem 1 Consider the Wigner ensemble of random real symmetric matrices with inde- 
pendent entries defined by flO|)-(tO|). Assume additionally that the third and fourth moments 
of Wjk do not depend on j and k and that JT^ = supj <k 'E{\Wjk\ 5 } < +00. 

Then the following asymptotic formula 



rrS}Hz) = r[z 




w 2 r 2 (z) ar A (z) 



[1 — w 2 r 2 (z)] 2 1 — w 2 r 2 (z) 
holds for any z G U (U is defined in ( |fl.22| ) ) 



+ 0(n' 3 / 2 ). (HI.31) 



Proof. By the resolvent identity fLT9|) - fL2QD , 



-z~ x + (zn)- 1 £ E{G jm H mj }. (111.32) 



j,m=l 



If we were following the conventional perturbational-diagrammatic approach trying to develop 
the asymptotic expansion for E{g n (z)}, we would repeatedly iterate the resolvent identity 
selecting on each step the terms that contribute to the leading and sub-leading terms. The 
obvious drawback of such approach is that infinitely many iterations are needed and in the 
non-Gaussian case, when there is no analogue of the Wick theorem, the diagrammatic approach 
is rather complicated. 

We propose making use of ( ll.l(j| ) instead of iterating the resolvent identity. For each pair 
(j, m), Gj m is a smooth function of H m j and its derivatives are bounded because of ( |II.18| ) and 
the inequality 

\G jm \ < \\G\\ < \Qz\~ 1 (111.33) 

which holds for the resolvent of any real symmetric matrix. In particular, \D^Gj m \ < C\Qz\~ 5 
where C is an absolute constant. Here and thereafter we use notation D m j for d/dH m j. 

According to (|T2|)-([L"3l) and our assumptions, the fifth absolute moment of H m j is of order 
n -5/2 ^hus applying (|II.16|) (with p = 3) to each of the summands in the r.h.s. of (|III.32|) one 
finds that 

3 1 n fl i S \(a+l)/2 

zmgHz) = -1 + E "7^3)72- E +l1 J E {D« mj G jm } + e„, (HI.34) 

a=l 71 j,m=l a - 

where the semi-invariants of W m j and 

C JT 5 



< 



n 3/2 |Q;^|5 

Obviously, G is a complex symmetric matrix, i. e. Gj m = G m j- By ([LL18Q , -D^ m G>, 
a!Giti and 



mm 



— ^injGjm = Gj m + GjjG mm (III. 35) 

D m jG jm = 2Gj m + 6Gj m GjjG mm (III. 36) 

D m jG j m 6Gj m -\- 3QGj m GjjG mm -\- QG jjG mm (III. 37) 

for distinct j and m. Let us set k<i = w 2 and k 4 = o in ( [111. 34] ) . Then, as a consequence of 

flnn^ )- flnnr?D , 

zmU(z) = -1 - W 2 m^, z) - n~ VERW} + aE{^(z)}] + e n , (111.38) 

where 

j n \ n 

c n (z) = - £ G 2 m , d n (z) = - E Cm (HI.39) 



n i,m=l n m=l 



and 



^3/2 (jC 

provided |S2|2w and n is large enough 



+ • ( IIL4 °) 



To infer ( [III.38|) - (|III.40| ) from ( |IH.34j) , notice first that for the sum over j = m in the r.h.s. 
of ( |Iii.34| ) we have the bound 



n (a+3)/2 



m=l 



< 



c 



1 



OC 



< 



1 



n 



(a+l)/2\Q z \a+l ~ n (a+l)/2 — n 3/2 



a = 2,3. 



for all realizations of Wjk- Therefore, being interested in the leading-order and l/n-order terms 
of (z) we can omit 5j m from the factor in front of the second and third derivatives. As for the 
first derivatives, it follows from ([III. 35| ), that for all j and m (l + 5j m )Dj m Gj m = G 2 m + GjjG mm 
and the term arising from 5j m contributes to 1/n-order term in the asymptotic expansion of 

Now, G 2 m in the r.h.s. of ( PX^ ) makes E{c n {z)} in (pTQgp and G^G mm does m^{z, z).The 
term containing G 2 jG 2 mm in the r.h.s. of (|III.37| ) leads to E{d 2 l {z)} in ( |III.38|) and the rest in 
the r.h.s. of and ( [111.37] ) contributes to e n in ( 111.38| ). Corresponding bounds for the 

terms coming from G^ m in ( |III.36|) and from G^ m and G 2 m GjjG mm in (|III.37|) result from the 
simple inequality 



j,m=l 



(111.41) 



which holds for the resolvent of any real symmetric matrix. Estimating the term coming from 
Gj m GjjG mm in the r.h.s. of ( |111.3(j| ) requires a longer calculation. Set 



h n (z) 



E GjmGjjG Ti 



j,m=l 



(111.42) 



Substitute the r.h.s. of ( [11.20|) for Gj m in h n (z). Then 

1 n 

zE{h n (z)} = — £ E{G 2 rnm \ - w 2 E{g n (z)h n (z)} + 0(n 



-l/2\ 



(111.43) 



m=l 



as it follows from ( |11.16| ), ( |111.35| )-( [111.37| ) and simple resolvent bounds like ( [111.33| ) or ( |111.41| ). 
Here and below we use notation 0(n~ p ) for remainders admitting the upper bound Cn~ p , where 
C does not depend on n for > 2w. 

According to (24) the variance of g n (z) is of order n~ 2 under our assumptions. In other 
words 

m£> (z, z) = [mf (z)] 2 + 0(n- 2 ) (111.44) 

if z G Uq. Obviously, 

E{h n (z)g n (z)} - E{h n (z)}E{g n (z)} = E{h n {z)[g n (z) - E{g n (z)})} 
and by the Cauchy-Schwarz inequality 

E{h n (z)g n (z)} = E{h n (z)}m^(z) + 0{n~ x ). 

Therefore by (|ITl~4l|) , 

[z - w 2 mW(z)]E{h n (z)} = -- ± E {G 2 mm } + 0{n^ 2 ) (111.45) 

n 771=1 

and E{h n (z)} is of order unity. The term we wish to estimate is 



E{Gj m GjjG mm } — — ^E{h n (z)}. 



n 



5/2 



j,m=l 



and from (|III.45|) we see it is of order n~ 3 ^ 2 . This proves (|III.38|) -( |III.40| ). 

The calculation above is typical of our approach and uses ( ITTB ) and ( |il.2U| ) combined with 
simple resolvent bounds on different stages. In what follows we shall ( |I14|) and ( |I.15|) often use 
similar calculations omitting details. 

Equations ( pL^) - (|ITL40|) and ^TT4|) imply that 



r(z) + 0{n 



(111.46) 



where r(z) solves the equation ( |11.30| ). Because of ( [T^ ) m l n (z) as a function of z must satisfy the 
inequality ^sr{z)^sz > 0. This restriction fix the branch of the square root in the expression for 
the solutions of (|H.30|) . Thus r(z) coincides with the Stieltjes transform flt.lOf) of the semi-circle 
law ( |1.5|) , as expected. 

Once the leading term of m$(z) is found, we can proceed with finding the sub-leading 
term. From (|IH.38| ) it is clear that performing this task requires calculating the leading-order 
terms of E{c n (z)} and E{d 2 l (z)}. Substitute the r.h.s. of (|II.20|) for one of Gj m in c n (z) and 
apply ( |11.16| ). As a result, 



zE{c n {z)} 



-m£\z) 



2w 2 E{c n {z)g n {z)} + Oin- 1 ). 



By (ITU), 

zE{c n (z)} = -m^\z) - 2w 2 m[ l \z)E{c n (z)} + O^ 1 ) 
and in the leading order 

E{c n (z)} = -mSHz)[z + 2w 2 m£\z)]-\ 
Taking into account ( |111.46| )-( |11.30| ) we conclude that 

E{c n (z)} = r 2 (z)[l - w 2 r 2 (z)}- 1 + O^ 1 ). 



(111.47) 



(111.48) 



Now, calculate the leading-order term of E{d 2 l {z)}. Recall that according to ( [II.24D the 
variance of g n {z) is of order n~ 2 . By (H.ll|) , g n {z) — n ~ X ELi G mm converges almost surely to 
r(z) as n — >• oo. Or, put it into another way, the Cesaro limit of G mm is r(z). This suggests 
that the Cesaro limit of G 2 nm should be equal to r 2 (z), or in other words d n (z) should converge 
almost surely to r 2 (z). Therefore E{d 2 l (z)} should converge to r 2 (z). 

To prove the convergence rigorously and to estimate its rate, we first note that the variance 
of d n (z) is of order n~ x if z G Uq (this can be proved following calculations of Appendix |Bj) . 
Therefore 



E{d 2 n {z)}=E{d n {z)} 2 + 0{n 



(111.49) 



Thus, it suffices to find the leading-order term of E{d n (z)}. 

Again, as in the case of c n (z), substitute the r.h.s. of ( |U.20p (j = m) for one of G mm in 
d n (z) and apply ( p,16| ). As a result, 



By ( ggp , 



zE{d n (z)} = -m£\z) - w 2 E{d n (z)g n (z)} + Oin' 1 ' 2 ) 



zE{d n {z)} = -m%\z) - w 2 mW(z)E{d n (z)} + 0{ n - l l 2 ) 



and 



(111.50) 



E{d n (z)} = -wg\z)[z + w 2 m^\z)\- 1 + 0{ n - 1 ' 2 ). 
Finally by (ITT061) and (gOg) , 

EK(z)}=r 2 (^) + 0(n- 1 / 2 ) 
and by jlL49|) , 

E{rf 2 ( 2 )} = r 4 (z) + 0(n- x / 2 ). (111.51) 

Now we are in a position to find the sub- leading term of m^'(z). Collect ( pi.38|) , ( |III.44| ), 
HLH) and (EIDfiD and write 



2m W( 2 ) = _l_ w 2r (l)/^|2 _ 1 

n 



w 2 r\z) 4 

+ err m 



1 — u? 2 r 2 (z) 



+ 0(n- 3 / 2 ) 



In view of ( |III.46|) and ( |II.30| ) this relation is obviously equivalent to the statement of the 
theorem, i.e. to the asymptotic formula ( [III. 31 ). The theorem is proved. 



Remarks 

1) Our bound for the remainder in ( |111.31[ ) is an optimal one. For assuming the 6th absolute 
moment of Wjk to be uniformly bounded and keeping one more term when applying ( |II.16 ), 



-3/2 



in the asymptotic expansion of m^'(z). This term is 



we can find the term of order n~ 
proportional to k 3 = E{W? fc } 

2) If the distributions of Wjk are such that E{PF 3 fc } = 0, then the bound 0(n~ 3 ^ 2 ) for 
the remainder in ( |H1.31| ) can be strengthened to 0(n~ 2 ). For terms of order n~ 3 ^ 2 appear in 
(|HI.31| ) due to contribution of K 3 n~ 3 ^ 2 E{h n (z)} to e n in (|III.38| ) and also because of ( |III.51| ). 
If E{W^} = 0, then k 3 = and we can prove that the remainder in (111.51) is of order rT l 
(terms of order n~ 1//2 in the r.h.s. of ( [III.51|) are proportional to Ks). 

3) For Gaussian Wjk the excess a is zero and ( |111.31| ) reduces to the asymptotic formula 



E{g n {z)}=r{z) 



1 + 



w 2 r 2 (z) 



+ 0{n~ 2 ). 



n [1 — w 2 r 2 (z)] 2 

that has been derived earlier by the formal diagrammatic approach!!. 



IV. LEADING ORDER OF F n ( Zl ,z 2 ) 

Let us recall our notation F n (zi,Z2) (see ( |L.14|) ) for the covariance function of g n (z) 



n-Hi (H-zI)- 1 , 

F n ( Zl ,z 2 ) = E{gC( Zl )gC(z 2 )} = E{g n { Zl )g^{z 2 )} 
In this section we prove the following 

Theorem 2 Consider the Wigner ensemble of random real symmetric matrices with inde- 
pendent entries defined by (^^)-^T3{). Assume additionally that the third and fourth moments 
of Wjk do not depend on j and k and that = sup J<fc E{|W^- fc | 5 } < +oo 



Let f(zi,z 2 ) be the function given by ( [i,15| ). If z\ and z 2 belong to U ( \ll.2'Jty , then the 
following asymptotic relation 



is valid. 



F n (z u z 2 )=n~ 2 f(z u z 2 ) + 0(n- 5 / 2 ) 



(IV.52) 



(IV.53) 



Proof. Let us first prove ( |IV.52|) under assumption 

Jpj = supE{\W jk \ 7 } < +00. 

Let Gj m (z) denote a matrix element of (H — zl)^ 1 . By (|11.20| ), 

1 n 

z 1 F n (z u z 2 ) = - J2 ^{H mj G jm (z 1 )g^(z 2 )}. 

71 j,m=l 

For each pair (j, m) Gj m (zi)g^(z 2 ) is a smooth function of H m j and its derivatives are bounded 
because of (|IlT33l) . In particular, \D 6 mj [Gj m {z 1 )g°{z 2 )\\ < C^z^ 1 + (|^ 2 | -1 ) 8 - Therefore 

by (grig), 



zF n (zi, z 2 ) = 

5 n 

X! r? (a+3)/2 X! 
a=l j,m=l 



(IV.54) 



Ka+l(l+5 im ) (a+1)/2 



E{C[G jm ( 2l K c ( 22 )]}+ £ni 



where K a are semi-invariants of W m j, as in (pil.34j ) , and 

kni^-^cvMi^ii^ + a^r 1 ) 8 

Performing differentiating in the r.h.s of ( |IV.54|) one finds that the sums of fifth, fourth 
and second derivatives in the r.h.s. of ( |1V.54| ) contribute to ziF n (zi, z 2 ) terms of order n~ 9//2 , 
n~'l 2 and n~ 5 ^ 2 , respectively (corresponding bounds can be obtained using ([11.24 ), ( III. 33 ) and 
(|IIL41| )). So, these derivatives give no contribution to the leading-order term of F n (zi,z 2 ). It 
remains to estimate the contributions coming from first and third derivatives. 

The contribution of third derivatives to ziF n (zi, z 2 ) consists of several terms which we shall 
label by integer a and h satisfying < a, h < 3 and a + b = 3. These terms are 

ka n 

s^ b \z llZ2 ) = ^ Y, V{D a m3 G 3m {zi)D b m3 g C n {z 2 )}. 



j,m=l 



First estimate s ( ^ ,0 \z\ 1 z 2 ). After differentiating it takes the form 

s^°\z 1 ,z 2 ) = -n~ 2 K i B{n- 1 ]T G 4 jm { Zl )g^(z 2 ) + 6h n (z 1 )g^{z 2 )} 

j,m=l 

-n- 1 K 4 E{d 2 n (z 1 )g^(z 2 )} (IV.55) 



(h n (z) and d n (z) are defined in ( |III.42|) and ( |III.39|) , respectively). As it follows from (|IL24|) 
and ([111.41 ), the mean value in the r.h.s. of the equation above is 0(n _1 ), so 



,(3,0), 



21,22, 



n-VEKOzOffe)} + 0(n 



Now we employ the obvious algebraic relation (in the below sub-index C indicates subtracted 
mean value) 



E{r, 2 ^} = 2E{r^'}EM + EWfe'} 



-PtC 



(IV.56) 



and write E{d 2 l (zi)g^ (z 2 )} as 



2ER(^(z 2 )}EK(z 1 )} + E{K(^)]^(^ 2 )}. 

If z G C/o, variances of g n (z) and o? n (^) are of order n~ 2 . In addition to this, d n (z) is bounded 
in absolute value by Cl^zl -2 for all realizations of Wjk- Therefore by the Cauchy-Schwarz 
inequality, E{<i 2 (zi)g%(z 2 )} = 0(n~ 2 ), provided Zi,z 2 G U . Thus we have proved that 
s^'°\zi, z 2 ) = 0(n~ 3 ). A similar argument shows that s£> 3 \zi, z 2 ) and s^ ,l \z\, z 2 ) are 0( 
too. The last term we need to estimate is s£' 2 >(zi, z 2 ). It is easy to see that 



n 



,(1.2)/ 



n 



n \ Z li Z 2 
-1 



2K A E{n 1 J2 G mm ( Zl )G mm (z 2 )}E{n 1 J2 G mm ( Zl )G 2 jm (z 2 )} + 0(n~ 3 ). 

j,m=l 



m=l 



Mean values in the above are calculated in exactly the same way as E{c n (z)} and E{d n (z)} 
have been done. For large n: 



n 



-^/^{rT 1 £ G mm (z 1 )G mm (z 2 )} = r(z!)r(z 2 ) + 0(n- 1 ' 2 ) 



m=l 



and 



E{n~ x £ G mm (z 1 )G J 2 m (^)}=r(^ 1 )r 2 (^)+0(n- 1 ) 

j,m=l 

(compare with ( |111.48| ) and ( |111.50| )). Thus we conclude that the contribution of third derivatives 
is 



1 err 2 \zi)r 3 \z 2 ) 
n 2 1 — w 2 r 2 (z 2 ) 

( |11-17| ) (we recall using o for « 4 ). 

First derivatives in the r.h.s. of (|IV.54j) contribute to Z\F n {z\, z 2 ) the term 

t n {zi,z 2 ) = -w 2 E{g 2 ( Zl )g%(z 2 )} - rr 1 w 2 E{c n (z 1 )g^{z 2 )} 
-n- 2 2w 2 E{n-Hr (H - zjy\H - z 2 iy 2 }. 

By the resolvent identity (HI. 19 ) 

(H - zJ)-\H - z 2 iy l = ( Zl - z 2 )- l [{H - zjy 1 -(H- z 2 iy l \ 

Thus one reduces E{n~Hr (H - ziiy l (H - z 2 I)~ 2 } to 

(zi - z 2 ) -1 [(zi - z 2 y l E{g n ( Zl ) - g n (z 2 )} - E{c n (z 2 )}]. 

Now recalling ( POSD , (gTrg ) and (|IL30|) , 

E{n~Hr (if - ztiy^H - z 2 iy 2 } = 
r(zx) - r(z 2 ) 12 



r(zi)[l — w 2 r 2 (z 2 )] 



Zl - z 2 



+ 0(n- 1 ). 



(IV.57) 



(IV.58) 



(IV.59) 



Clearly, E{c n (zi)g%(z 2 )} = E{c^(zi)g^(z 2 )} and the corresponding summand in the r.h.s. 
of ( [IV. 5^ ) is 0(nr 3 ). So it remains to find E{g 2 (zi)g%(z 2 )}. 
Use (|iV.56|) to write 



E{^i)<£(* 2 )} = 2F n ( 2l)22 ) m W(, 1 ) + E{[^(, 1 )]\ C (. 2 )} 
= 2F n (z 1 ,z 2 )m£\z 1 ) + 0(n- 5 / 2 ) 



(IV.60) 



The latter equality uses E {^n(^i)] 2 ^n(^)} = 0(n" 5/2 ), the bound which can be obtained 
following calculations of Appendix |B|. 

Now we are in a position to find the leading order of F n (zx, z 2 ). Collecting ( PV.57|) -( |TV60|) , 
we find that 



ZiF n (z 1 ,z 2 ) 



-2w 2 F n (z l ,z 2 )m^(z 1 



2w 2 



1 ar 2 (z 1 )r 3 (z 2 ) 

n 2 1 — w 2 r 2 (z 2 ) 
2 



r(zi) - r(z 2 ) 
z\ - z 2 



r(^i)[l — w 2 r 2 (z 2 )] 
As it is clear from ( PTL46D and ( |1L30| ), 

[ Zl + 2w 2 m { ^\z l )}- 1 = -r(zi)[l - w 2 r 2 {z l )}- 1 



+ 0(n" 5 / 2 ) 



0(n 



and we end up with ( PV.52|) . 

The standard truncation technique of probability theory allows to prove (|IV.52| ) in the case 
when only 5th absolute moment of the random variables Wjk is uniformly bounded. Calcu- 
lations using the truncation technique are similar to those used in next section in proof of 
Theorem 3 and we omit them. Theorem 2 proved. 

One can consider the covariance function F n (zi,z 2 ) for the Wigner ensemble of random 
Hermitian matrices (see remark 3 after the statement theorem 3 in next section). Repeating 
almost literally calculations used in proof of theorem 2, one can prove that for the Wigner 
ensemble of Hermitian matrices ( |lV.52j ) is still valid. The only difference is that now f(zi,z 2 ) 
is given by the r.h.s. of Eq. (|I.15| ) multiplied by factor 1/2. 



V. GAUSSIAN FLUCTUATIONS OF THE CENTRALISED TRACE OF THE 

RESOLVENT 

In this section we prove the statement which is analogous to the central limit theorem in 
the same sense in which the result ( |I.11|) is analogous to the law of large numbers. Indeed, we 



can rewrite ( |1.11| ) as following limiting relation 

1 n 

lim - E G mm = r(z) (V.61) 

valid with probability 1. Since the l.h.s. here has the form of the arithmetic (Cesaro) mean, 
this relation is obviously similar to the strong law of large numbers (or more generally to the 
ergodic theorem). Common wisdom of probability and ergodic theory suggests that (|V.61|) 
should imply that the probability distribution of the random variable 



n 1 ' 2 



1 E ( G ™ m - E{G mm )}l = n 1 ' 2 ^) - E{g n (z)}} (V.62) 



n 

771=1 



has the Gaussian form in the limit n = oo. We prove that under rather natural conditions on 
Wjk this is indeed the case provided that we use non-standard normalisation, replacing n 1//2 in 
( |V.62|) by n, i. e. we consider just the centralised trace of the resolvent 

l in \z) = E ( G mm ~ E{G mm }) = ng n (z) - E{ng n (z)} (V.63) 

m=l 



instead of n l ^ 2 ^ n \z). This normalisation can of course be anticipated from the formula ( [1.14] ) 
giving the order of magnitude (in fact, the asymptotics) of the variance of g n {z). This decay of 
the variance, which is "twice" more strong than in the standard central limit theorem setting, 
is rather typical for a number of problems of the theory of disordered systems with non-local 



interaction and is known as the strong self-averaging property (see e. g. Refs. 28 and 33 



Theorem 3 Consider the Wigner ensemble of random real symmetric matrices with inde- 
pendent entries defined by flL^)-(|0|) assuming additionally that the fourth moments of Wjk 
exist and are independent of j and k and that the probability distribution functions Pjk(w) of 
Wjk satisfy the condition: for any fixed fixed v > 

Jim — J x 2 dProb[W jk < x] = 0, for any v > 0. (V.64) 

j^k II 



Then for any z from Uq = {z G C± : \0>z\ > 2w} the random function ^ n \z) (\V.63j ) converges 



in distribution as n — > oo to the Gaussian random function 7(2) with zero mean and the 
covariance function f{z\,Z2) given by (jTl5|). In other words, for any integer q and arbitrary 



collection z±, . . . , z q of complex numbers from Uq the joint probability distribution of random 
variables 7^(21), . . . ,^ n \z q ) converges as n — > 00 to the q-dimensional Gaussian distribution 
with zero mean and the covariance matrix [f(z s , z t )] q st= i 



Remarks 



1. Limit theorems concerning ^ n \z) for the Wigner ensemble were established for the first 



time by Girko (see Ref. [24] and references therein) under assumption that there exist a positive 
S such that 



sup E\W jk \ i+d < 00, (V.65) 

which is slightly more restrictive than (|V.64|) . For example in the case of identically distributed 
Wjk, ([V-64 ) is obviously satisfied if W4 = E-fWA} is finite. However, the more important in 



our opinion improvement of the result of Ref. |23 is that we calculate the covariance matrix of 



the limiting Gaussian process in the explicit form while in Ref. 24 this matrix was given in the 



implicit form as a solution of a system of cumbersome partial differential equations. 

2. For the random variables Wjk satisfying (|V.65|) we can estimate the rate of convergence: 

sup I E{^ n \ Zl )^ n \z 2 )} - f{ Zl ,z 2 )\ =0(n^ 2 ). (V.66) 

zi,z 2 €Uo 

3. Consider the Wigner ensemble of the n x n random Hermitian matrices defined as in 
(|L|) with Wjk = Ajk + iBjk, j < k, Wjk = WL, where Ajk and Bjk are mutually independent 



random variables with zero mean, variance w 2 /2 and excess a/2. It can be proved by analogous 
technique that in this ensemble the fluctuations of the trace of the resolvent around its mean 
become Gaussian in the limit n —>■ 00. The corresponding covariance function is given by (|L15|) 
in which the factor 2 is replaced by 1 in the denominator of both terms. 



Proof. We shall work with real- valued variables a^ n \z) = $l^ n \z) and f3^(z) = ^sj^(z). 
Then we have to prove that the limiting random functions a(z) and (3(z) are jointly Gaussian, 
i.e. if 



X(z,c) 

and 

(a(c),b(c)) = 



a(z) if c = a; 
(3(z) if c = 0, 



{ (1/2, 1/2) if c = a 
\ (l/2z,l/2z) if c = /3 



then E{X(,2,c)} = and for any integer g and arbitrary collections z s ,s = l,...,q, z s G Uq 
and c s , , s = l,...,q, c s G {a,/3} the joint probability distribution of random variables 
X(z 1 ,ci), . . . ,X(z q ,Cg) is the q- dimensional Caussian distribution with zero mean and covari- 
ance matrix 

E {X(z s , c s )X(z t , c t )} = a(c s )a(ct)f(z s , z t ) + a(c s )b(c t )f(z s , z\) + 

a(c t )b(c s )f(4, zt) + b(c s )b(c t )f(4, 4), (V.67) 

Let us consider the characteristic function of random variables X(zi,Ci), ..,X(z q ,c q ) which we 
shall write in the form 



Q 

E 

8=1 



e q n \T g ,C q ,Z q ) = E{nexp{zr s [a(c s ) 7 (n) (^) + Kc,)7 (n) (4)]}, 



where T q = (n, . . . ,r q ), C q = (c x , . . . , c q ), Z q = (zi, . . . , z q ) 

Recall that we designate the complex conjugate by the symbol f. Also writing the char- 
acteristic function we shall often omit indices indicating its dependence on n and some other 
variables provided there will arise no confusion. 

Obviously 

^ E{e q (T q )} = iB{e q [a(c s )^(z s ) + b(c s )^\zl)}}. 

Our aim is to show that there exists a set of the "covariance" coefficients , s,t = l,...,q 
such that for each fixed T q 



lim 



E{e<>W«)(^) + 6 s7 (n) (4)]} -^4? E{e^} 

t=i 



0, z G U , 



to show that limits of all these coefficients exist 

A st = lim 4?, (V.68) 

n— >oo 

and correspond to the r.h.s. of ( |V.67| ). Then standard arguments will allow us to prove that 
the limit characteristic function has the Gaussian form exp(— 1/2 J2t t=i A s tT s Tt). 
Thus, we have to compute 

n 

E{e ?7 (n) (*)} = £ K{e°{Z q )G 3j }, 
j=i 

for large n (we recall that sub-index C indicate centering to zero mean). Then putting one of 
Zi, . . . , z q or one of their conjugates in place of z we calculate the limits in ( |V.68| ) 

We have complex energies z, Zi, . . . , z q and we introduce notation G(z s ) for the resolvent 
corresponding to z s keeping notation G for the resolvent corresponding to z. 

By the resolvent identity ( |11.2C| ), 



E E{e^,} 



z 1 E E{e^G im if mi }. 

j,m=l 



(V.69) 



We compute the average in the r.h.s of ( |V.69|) following the scheme described in Section 2. 
However its direct application requires too strong conditions on the distribution of Wjk- Thus, 
we modify slightly the general scheme and carry out more accurate estimates. 

Denote by E mj - the conditional expectation Ej-lPVmj = w} and rewrite the right-hand 
side of ( |V.69| ) in the form 



n „ 

J E mj {e^G jm }wdP mj (w). 



j,m=l 



We split the integral into the two ones over sets Ti(n) = {u : \W m j\ < 5n 1 / 2 } and T 2 (n) = {u : 
\W m j\ > 5n 1 / 2 }. Now the inequalities 



n 



-1/2 



I E mj {e^G jm }wP mj (w) <r)- x n 1/2 £ / \w\dP rnj (w) 

n „ 

< r'^n- 2 / w A dP mj {w 

1 JToM 



j,m=l 



and the assumption ( |V.64j ) imply that only the integrals over Ti(n) gives the a non- vanishing 
contribution contribution to ( |V.69[ ). 

Following our general scheme, we expand the function e^Gj m in powers of random variable 
H m j — n ~ l ^ 2 W m j restricted to Yi[n). Since it is bounded in absolute value by v we can write 
the relation 



n 



-1/2 



n „ 5 

Y / ~E m j{eqG jm }wdP mj (w) = E S, 

j,m=l JFl k=l 



kin), 



(V.70) 



where 



S k (n) = n k ' 2 Y E 



e c G- 



j,m=l 



(*-l) 
mj 



w dP mj (w), 



and denotes that the k-th derivatives with respect to H m j is taken and then H m j is 

replaced by zero. Let us note also that in 65 the expression in square brackets is taken at some 
point H mj e (0, v). 

The term S\{n) vanishes as n — > 00 due our assumption (|V.64|) : 



\Si(n)\ < 



n 



-1/2 



E e 



j,m=l 



e c G- 



wdP mj {w) 



< rj V 3 n 2 



Y / w A dP mj (w) 



j,m=l JTl 

where r\ = \$sz\. 

The term S5 (n) vanishes as n — > 00 because it can be estimated by 



B 4 (T P 



E 



r? 6 n 5 / 2 . . 



w\ 5 dP mj (w) < 



B 4 (T p )v 
nfn 2 



Y / w 4 dP mj (w) 



where B±(T p ) is the upper bound of absolute value of the fourth derivative in (|V.70|) . For any- 
fixed T p , B^iTp) is finite and recalling ( |V.64[ ) we see that S§{n) goes to zero as n — > oo 

The term S3 (n) also vanishes as n — > 00. We establish this fact at the end of the proof. 

Terms S^n) and S^n) give main contribution to (|V.69|) . Let us first consider S2 (n). The 
resolvents G(z s ) and G are complex symmetric matrices and we have: 



E $G jn 

j,m=l 
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P c r r 

q M mm 



rnj 



n 
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j',m=l 
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n 

< E 

j,m=l 



y 

^ 2r s (a s [G 2 (> s )] jm + 6 s [G 2 (4)]jm) G jm e G 



.5=1 



(V.71) 



Each term of the right-hand side of this relation is a function of H in which if mj is replaced 
by zero and we have to come "back" to expressions dependent on the whole matrix H . To this 
end, we use again the resolvent identity but now in the "opposite " direction. We obtain for 
the first term of ( |VTTD 



w 2 n 1 



E E{ 

j,m=l 



q ^- r jj^ J mm 



>>'.) 



} = v * n - 1 ]T E { e qGjjG mm } - w 2 ¥ n \ 

j,m=l 



where 



!jr(n) 
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E E{ 



j,m=l 
n 

2^n 2 £ E{ 

j,m=l 



( } / wdP mj (w) 



Cq Cjj G mm 



( ) } / w 2 dP mj {w) 



i,m=l 



Co Gmm 



(3) 



W j. 



It is easy to see that the first and the last terms of \l/^ n - ) vanish as n — > 00 due to our assumption 
Using ( |V.64j ) and (|IH.37| ), we can rewrite the second term of in the form 



n 



j,m=l 



q 

e g E r * (a s [G 2 (z s )] mm [G(z s )]j, + b s [G 2 (z ] s )] mm {G(z\)\ j ^ G j:j G r , 

s=l 



where the remainder $( n ) includes terms which have one or more factors Gj m or terms of the 
form 



n- 2 E E{[e^[G(^y 2 [G( 

j,m=l 



It is clear that in all these expressions we can remove square brackets [■■■} m j because this 
procedure will add terms of order 0(n -1 / 2 ) to the sums under consideration. Using the estimate 
( |III.41| ), taking into account the self- averaging property 



E \g^(z)\ 2 = o(n 1 ), as n — > 00 
(see Lemma 1 of Appendix [B] for the proof), and relation 



lim E I (n- 1 £ G£ G^n" 1 £ G mm G 



' mm) 



(V.72) 



(V.73) 



with some a, (3, fi, v = 0,1,2 (Lemma 2 of Appendix |BJ), it is easy to prove that $( n ) also 
vanishes as n — > oo. 

We obtain finally that among terms coming from first summand in the right-hand side of 
( |V.71| ) only the following 



(V.74) 



— 2iwi 



E {ej £ r s . E { £ (a s [G 2 (^)] mm [G(^ s )] ij + 6 s [G 2 (z s t )] mm [G(zi)]„) G„G mm }. 
i,fc=i 



3=1 



does not vanish as n — > oo. 

The second summand in the r.h.s. of (|V.71|) vanishes as n — > oo. This becomes clear 
after applying the same procedure of removing square brackets [. . .] mj - to the expression 



n 



v2l(*) 



" 2 E™, m =i[ e f ( G im) 2 ]mj and usm S the estimates (|111.41|) and Q7773D . 



Let us consider the contribution of the last summand in r.h.s of (|V.71|) for a fixed parameter 
z s , s — 1, q. Taking into account ( |V.64| ) and repeating the "returning" procedure, we obtain 
for this term 



^ \[G ( z s)]jmG m je q 
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.2) / w 2 dP mj (w)-n- 5 / 2 E { [G 2 (^)] J - m G , mi e ? "W^. 

Fl i,m=l 



It is easy to see that in this equality terms with the first and the third derivatives vanish 
as n — > oo. The second derivative gives 



{\G (Zs)]jj Gjj G mm G mm -\- [G (z s )] mm G mm \G(z s ^jj Gjj^&q 



illj 



and 24 terms having a factor of the form (G a )j m , a = 1, 2. Omitting brackets [■■■] m j an d using 
( |V.72|) -( |V.73|) , we see that the last term of (|V.71|) gives the leading contribution 
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(V.75) 
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Consider now the term 64(71) of (|V.70|) . The third derivative [e^Gjm]^- consists of 140 
terms. One part of them vanishes as n — > 00 due to the property ( |V.72| )-( ]V.73D , another part 
- due to the presence of the factor of the form (G a )j m , a = 1, 2, 3. Only six terms of the form 
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are non- vanishing in the limit n — > 00. These terms arise when we differentiate Gj m once and 
e q twice with respect to H m j. Combining these terms with ( |V.74j ) and ( |V.75| ), we finally obtain 
that 



E {e^tr G} = % — 1 E {e g c } £ t.^tT 1 E {tr (a s G(^) 2 + h s G{z\f) G} + (V.76) 
z zw g n yz) 
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mm b s G [z g ) jj G(z^ mrn jGjjG mrn }) . 

j.m=l 

Notice, that the denominator in the first term of the r.h.s. of this expression is bounded 
away from zero because z belongs to the domain (|II.22|) . Now, combining ( |I.11|) and ( |I.10|) with 
relations 

m 

Jim E^D^W]^} = r 2 (z s )[l - ^ 2 r 2 (z s )]-V(z), (V.77) 
j'=i 

we derive from ( |V.76| ) the final form of the covariance. 

Relation ( |V.77| ) can be easily deduced from our proof of (|V.72[) - ([V.73[) . 

Let us briefly discuss now the proof of the fact that S 3 (n) of ( |V.70| ) vanishes as 



n 



oo. The second derivative [e^G-,™,]^ gives terms each having the factor of the form 
[(G a )j :/ (G !fc )j m (G /3 ) mm ] m j. The brackets can be simply omitted because the "returning" proce- 
dure adds terms of order 0(n~ l l 2 ). Now, regarding n~ l l 2 (G a )jj as vectors and (G m ) m j as the 
kernel, we can write inequality 

| n - 3 /2 j2 (G a ) ]3 (G m ) m] (G?) mm \ < r «-P- m n- 1 ' 2 

j,m=l 

which completes the proof of Theorem 3. 



VI. SCALING LIMIT AND UNIVERSALITY CONJECTURE 

We have presented above the rigorous derivation of asymptotic corrections (in fact expan- 
sions) for moments and more complex quantities constructed from the traces of the Green 
functions of the Wigner random matrix ensembles. Now we use our result to draw certain non- 
rigorous conclusions on the form of the leading term of the correlation function S n (Ei, E 2 ) of the 
formal level density p n {E) = n _1 tr 5(H — EI). Since p n {E) = N' (E) where N n (E) is defined 
in ( |1.4| ), then basing on the relation ( |1.7| ) one can conclude that the number of eigenvalues lying 
inside the interval (E 1} E 2 ) with the center at E will be N(E 2 ) — N{E\) ~ np n {E) ( y E 2 — Ex), i.e. 
that the mean distance between levels is [np n (£')] _1 . Thus the scaling E 2 — E\ = 0{n^ 1 ) defines 
the microscopic or local regime in which one deals with a finite numbers of eigenvaluesEI'Q. 

Consider the density-density correlation function 

S n (E h E 2 ) = E {pniEMEh)} - E K(£i)} E {p n (E 2 )}. (VI.78) 
By using (8) and (9) we obtain from (14) and ( |V.63[) that the Stieltjes transform of S n (Ei, E 2 ) 

F n ( Zl , **)=[[ Sn{EiE 2 ) dEidE Q ^ Q 
J J {Ex- z l ){E 2 - z 2 ) 



is 



F n (z 1 ,z 2 ) = n- 2 E{ 7 W(zi) 7 W(^)}. 

It follows from the inversion formula for the Stieltjes transform f(z) — j(E — z)~ l p(E)dE 

p(E) = 7r~ 1 ]im^f(E + ie)=I El {f(z)} (VI.79) 
40 



that to find S n (Ei, E 2 ) , one has to know F n {zx, z 2 ) up to the real axis in both variables because 



S n {E 1 ,E 2 )=I El cI E2 {F n {z 1 ,z 2 )} (VI.80) 



On the other hand, we have found the form ( 1.14 ) and ( |1.15| ) of F n (zx, z 2 ) only in the domain 



\$$z\ > 2w. However, since the function f(z\, z 2 ) given by (|[.15|) can obviously be continued up 
to the real axis with respect to the both variables Z\ and z 2 we can apply to the first term of 
(|L14|) the operation Ie x Ie 2 i E\ 7^ E 2 to compute formally the "leading" term of the density- 
density correlation function. This means that we perform first the limit n — > 00 and then the 
limits ei, e 2 | 0. This order of limiting transitions is inverse with respect to that prescribed by 
the definition of this correlation function. 

To make these computations , we use the identity 

r 1 - r 2 rir 2 



Z\ — z 2 1 — w 2 r 1 r 2 

which follows from ( |I.10|) or (|II.30|) . The identity yields the relations e\r(E + ie)\ 2 = ^sr(E + 
ie)(l — w 2 \r(E + ie)\ 2 ) and \r(E + i 0)| 2 = w~ 2 for E such that 5sr(E + iO) > 0. Combining 



these relations with (|L5|) , we obtain that 



w 2 [^r{E + iO)} 2 = — - and w 2 [Qr{E + iO)] 



Aw 2 Aw 2 



Using these equalities, we derive from our result ( [L.14 ) and ( |L.15|) and from ( |V1.79| ) and 
(|VL80D that 



/3ir 2 [n(Ei - E 2 )} 2 (Aw 2 - Eiyi 2 {Aw 2 - E 2 ) 1 / 2 

a (2w 2 - E 2 )(2w 2 - E 2 ) 

2n 2 7r 2 w 8 {Aw 2 - E 2 y/ 2 {Aw 2 - E 2 ) 1 / 2 ' 



(VI.81) 



with (3 = 1. It can be shown that for the Hermitian matrices with independent entries (see 
Remark 3 to Theorem 3) the density- density correlator has the same form with (3 = 2. For the 
Gaussian orthogonal and unitary ensembles (GOE and GUE) a = 0, and we recover the result 



P7i 2 [n(E 1 - E 2 )] 2 {Aw 2 - E 2 y/ 2 (Aw 2 - El) 1 / 2 



obtained in Ref. |33] and Ref. |35| 

We see that in a general non-Gaussian case the respective expression depends not only on 
the second moment of entries , but also on their fourth moment via the excess a. 

The remarkable fact is that this dependence vanishes in the microscopic (called also scaling) 
limit 

E h E 2 -> E, n(E 2 -Ex) (VI.82) 
Indeed, it easy to see that in this limit we obtain from ( |VI.81|) very simple expression: 



t lim S n (Ex,E 2 ) = --^—. (VI.83) 

n{E 2 -E 1 )^s (3-K 2 S 2 

According to Wigner and Dyson (see e.g. Ref. |l|), the exact large-s asymptotics for the limiting 
correlation function of the Gaussian ensembles are: —1 / (vr 2 s 2 ) (GOE) and — sin 2 irp(E)s/ (ti 2 s 2 ) 
(GUE). Comparing these expressions with our results, we see that the procedure of computing 



of the correlation function yields for the general case the expression coinciding with the large- 
s asymptotics of the Gaussian ensembles correlation function smoothed over energy intervals 
whose length is much smaller than the macroscopic scale w = E {W /2 } 1//2 but much bigger than 
the microscopic scale given by the mean level spacing [np(E)]~ 1 . It is natural to think that 
in our computations the smoothing has been implemented "automatically" due to the nonzero 
imaginary part of the spectral parameter Qzj. We notice that the same procedure is widely 
used in the mesocsopic calculations based on the Kubo formula, weak disorder perturbation 
theory, etc. 

The independence of the scaling limit expressions ( |V1.83| ) on the excess a can be regarded 
as a support of the universality conjecture for the Wigner ensembles. Let us mention supports 
of this conjecture for other ensembles. 

The first on concerns the so-called sparse (or diluted) random matrices whose entries are 
independently distributed random variables such that Pr {H^ i = 0} = p/n. The authors of 
Ref. [18| used the Grassman integral technique and found the Wigner-Dyson universal form of 
the density-density correlator if p is large enough. 

The secondly concerns the ensemble H = Y?u=i t m('' £ M )£ M 5 where r M and £ M = . . . , 
are independent identically distributed random variables (the ensemble was introduced in Ref. 
fZ7|). For this ensemble, whose entries are dependent random variables, the analogue of ( |V1.81 ) 
is obtained and it is shown that its scaling limit is the same as above. 

The third follows from Appendix |A] below. We consider there case of the deformed GOE 
(see definitions below). For this ensemble the analog of Theorem 3 was proved in Ref. In 
the appendix we present a short derivation of the density-density correlator and show that in 
the scaling limit it has the form ( |V1.83| ). 

We mention also that for the unitary invariant ensembles of the form (|L~T|) the universality 



conjecture is rigorously proved in Ref. [16] for a rather broad class of functions F(x) 
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APPENDIX A: SCALING LIMIT FOR THE DEFORMED GOE 

In this Appendix we find the exact form of the leading term of the covariance function 
F n (zi, z%) ( |L.14j ) for the ensemble Hd = + H, where if < -°- ) are n x n nonrandom matrices 
such that there exists the "unperturbed" IDS 

N^(E) = lim ra _1 #|e,- : e,- is an eigenvalue of and e,- < E} 

n— >oo 

and H belongs to the GOE ((|1.1|) with F = x 2 /Aw 2 ). This ensemble is calledi the deformed 
GOE. Because of the orthogonal invariance of the GOE distribution, we can restrict our con- 
siderations to the case of diagonal H^°>. So we assume that = [^fcej]™ fc=1 and real numbers 
Cj are such that the limit 

i n i 

g<®(z) = lim - Y —— (A.l) 



exists for all non-real z. The function g^\z) is the Stieltjes transform of N^(E). We shall 
use notation dj(z) for (e 3 - — z)~ l and gf n (^) for the normalized trace of the resolvent of H^. 
Subsequent arguments are quite similar to those used in derivation of ([II.24|) and (|II.26|) - 
.29|) for the GOE. By using (17) and (19) for Hi = and H 2 = one can derive the 



(A.2) 



following two relations ( analogues of (|II.26|) and ( 1I.27| ) ) : 

E{G jk (z)} = d j (z)5 jk + E{g n (z)G jk (z)}d k (z) + n~ 1 J2' E {Gj m (z)G km (z)}d k (z) 



and 



E{g c J( Zl )G n (z 2 )} = u; 2 E{^(^)G„(z 2 )}E{^(z 2 )}rf,(z 2 ) + 
w 2 E{g^z 1 )g^(z 2 )G jj (z 2 )}d j (z 2 ) + 
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'■n- 1 E{g^z l )[G 2 (z 2 )} JJ }d J (z 2 ) + 



2w 2 n- 2 Y,V{iG 2 (zi)} jm G m] (z 2 )}d J (z 2 ) 

m 

It follows from ( |A.2| ) that if for z e Uq where Uq is defined in ( [ll.22j ) 

lim E{\g n (z)-E{g n (z)}\} = 0, 



(A.3) 



(A.4) 



then lim^oo E{g n (z)} = g(z), where g(z) is the unique solution of the functional equation^ 



g(z)=g^(z + w 2 g(z)) 

satisfying $$g(z)$$z > 0. In the equation above, g^°\z) is given by ( |A.1| ). 
It is easy to show that ( |A.2j ) and ( A.4 ) imply the relation 

» 



(A.5) 



sup \E{G jj (z)}- 9 p(z)\ = 0(n- 1 ), z E U , 



j=\,...,n 



(A.6) 



where g\ n \z) solves the equations 



9i ( z 



ej — z — w 2 g^(z) 



, j = l,...,n gW(z) = n- 1 £ g 



(»)/ 



(A.7) 



m=l 



Indeed, if Vf> = E{g^(zx)G jj (z 2 )} then by (pOj ), 
vj n) = w 2 v} n) E{g n (z 2 )d j (z 2 )} + 



n 



1 J2 V^E{g n {z 2 )d 3 {z 2 )} + 2w 2 n' 2 E{J2[G 2 (z 2 )} :im G m] (z 2 )d 3 (z 2 )} + 



(A. 



Now, repeating arguments used at the end of Section II, we can easily obtain the estimate 
n ~ l J2j Vj — 0(n~ 2 ) which proves (|A.6|) . Using this estimate and considering ( A.8 ) once 
more, we obtain the estimate 



sup | Vj | 



0{n" 2 ). 



(A.9) 



It follows from the resolvent identity (|I1.19| ) that 

_ [G 2 { Zl )\ j3 [G{ Zl )} 3] - [G{z 2 )} 



Y,[G 2 (z l )} rm [G(z 2 )} 



jj 



rnj 



Zl ~ Z 2 



{zi - z 2 y 



Taking into account this relation, (|A.6|) , and (|A.9|) , we obtain that if 



then 



f d (z u z 2 ) = lim n 2 E{g c \z 1 )g u (z 2 )} 



\n 



- 1 J2V^-f d (z 1 ,z 2 )\ = o(n- 2 ), zeU , 



and as the result 

fd(zi,z 2 ) = 



2w2 lim - 1 ^^i) _ 9tK^)W{z 2 ) 



{Z X - Z 2 f n-Xio j _ w 2 n -l ^. g${z 2 y 



lim 

n— >oo 



n- 1 E i [G 2 (^ 1 )](^)^)(z 2 ) 



Zl - Z2 n^oo (1 _ ^-1 E . 0&)( Zl )2)(l _ W 2 n -1 £ . ^(^)2) 



(A.10) 



Since lim^oo n _1 X} m 9m\ z ) 2 = I(E — z — w 2 g (z))~ 2 dN (E) = $ 2 , then the second fraction 
of the last term of ( |A.10| ) is not singular for Zi = E ± iO with E such that ^sg(E + iO) > 0. 
Thus, this term vanishes in the scaling limit ( |V1.82 ). 

Consider now the first term of the right-hand side of ( |A.10|) . Simple computation shows 
that 



n 



-i 



E^ n) (^i)£ n) (^) 



g^(z 1 )-g( n \z 2 ) 



z x -z 2 + w 2 \gW ( Zl )-gW {z 2 )\ 
Since, according to (|A.4j) - ( |A.7| ) lim^oo (z) = g(z) we find for this term 
2^ 2 / 9a(zi) - ga(z 2 ) 



[zi - z 2 ) 2 \ zi - z 2 + w 2 [# (£i) - 9o(z 2 )} 



$ 2 (l-^j^ + o^-zai- 1 ; 



This relation implies that in the scaling limit QVI.82|) we obtain again the simple universal 
expression (|VI.83| ) . 



APPENDIX B: AUXILIARY FACTS 



Lemma 1. Self- aver 'ageness property (\V-7Q) holds under assumptions of Theorem 3 
Proof. We denote 

F n (z,z') = E{g°(g') c } = E{g c g'} (B.l) 

where g = g n {z) and g' = g n {z'), g n {z) = n~ x tr (H - ziy 1 and g c = g - E{g}. 
Obviously, G dj (zt) = GUz) and F n (z,z*) = E{\g c (z)\ 2 } . 



Let us apply the resolvent identity ( |H.20| ) to the last factor g' in the right-hand side of ( |B.1| ). 
We obtain the relation 



F n (z,z') = —J2E{g c G' jm H mj }, 



z'n 



(B.2) 



where G' = G(z'). Comparing relations (p.2j ) and ( |V.69| ), we see that their right-hand sides 
are similar. The only difference is that the sum in ( |B.2[ ) has extra factor n" 1 and e q of ( |V.69[ ) is 
replaced by g. Hence, one can compute the average in ([B.2|) in the same way as it was done for 



the right-hand side of (|V.69|) and come to the expression F n (z,z') = J2l=iTk(n) where Tk(n) 
are similar to Sk(n), k = 1, ..,5 in ( |V.70| ). 

Thus we find that T\{n) and Ts(n) are of order o(n _1 ), as n — > oo just as in the case of 
Si(n) and S3 (n). 

Consider T 5 (ri) which is analogous to £5(71) in ( |V.70| ). It contains four derivatives of GkkGj m 
by H m j. It follows from ( [III.37|) that the result of differentiating includes at least one factor 
Gj m . Combining ( |111.41| ) with inequalities used to estimate S 5 (n), we easily derive that T 5 (n) 
is a quantity of order 0(n~ 3 ^ 2 ). 

Let us estimate T±{n) acting in the same way as in the case of S±{n). As it was mentioned in 
0, the no n- vanishing contribution to S±{n) comes from terms arising from one derivative of G'j m 
and two derivatives of e q . The rest of the terms are of order o(l). Thus, in the corresponding 
terms of T 4 (n) we have to take into account only terms with factors G or G 1 having coincident 
arguments. It is easy to see that due to extra factor n _1 in front of the whole sum and factor 
rT l in g(z), these terms are of order n~ 2 . Thus, T^{n) is of order o(n~ l ). 

Turning to T 2 {ri) and taking into account previous arguments, we arrive at the relation 



FJz, z) 



^ £ E{g c G' j3 G' mm } - £ E{g c G' jm G' jm } + $'(,, z') 



(B-3) 



j,m j,m 

where $>'(z, z') = o(n _1 ). Using (|III.41|) , we easily obtain that 
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n- 



E i9 C ' G 'jm G 'jm} 



3,m 



<^EE{ij c ii:iG; m i 2 } 

j m 

"^ 1/2 {\9 C \ 2 } 



< 



n\^sz'\ 2 



Observing that 

E{g C g'g'} = 2E{g c g'}E{g'} + E{g c [g'\ c g'}, 
we derive from ( |B.3|) that for z' = z' and z eUq ( |II.22| ): 

C 1 F n (z,z^) - C 2 n- 1 \F n (z,z^)\ 1 / 2 - |$' n | < O, 



where C\ and C 2 are absolute constants (cf. (p.28|) . This inequality implies ( |V.72|) . Lemma 
proved 



Lemma 2. Under assumptions of Theorem 3 the relation (\V.73j is true. 

Proof. It suffices to show that 

R n = E{\G c 2 \ 2 } = o{\), n^oo, 



(B.4) 



where G 2 = n 1 J2jG 2 j. Repeating computations of previous proof, we obtain the following 



relation 



R n — — ^E{G 2 G' jj G' jm H rnj }. 



z n 



3,m 



Comparing again the right-hand side of this equality with those of ( |B.2| ) and ( |V.69| ) and re- 
peating the corresponding computation, we conclude that 



Rn ~ 'l/n 2 ^ E ^ G ^ G 'jJ G 'jj G 'mm} + H E { G 2 G 'jj G 'jm G 'mm} + 



3,m 



w 

z'n* 



(B.5) 



where $" n = o(l) as n — > oo for 2;' = ^ an( } > 0. Taking into account ( |V.72| ), we derive 
from ( |BT5| ) that for z e U 

C 3 R n - C 3 Rn - n- l C±\R n \ l/2 + o(l) < 0, 

where C3 and C4 are some absolute constants (cf. QlI.28|) ). This inequality implies ( |B.4|) . Lemma 
is proved. 
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